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1. Introduction 

The JLQCD collaboration has started a large scale simulation of dynamical overlap fermions [[I]], 
aiming at studying the QCD dynamics in the presence of very light quarks. For the overlap fermion 
formulation, large computational costs and non-locality of the overlap-Dirac operator have been 
practical and theoretical concerns. The computational cost is closely related to the near-zero mode 
density of the Hermitian Wilson-Dirac operator, 

H^{-mo) = rsDy^i-ma) = Y5{-mo + D^{0)) (mo = l+5), (1.1) 

where D^{0) is the standard Wilson-Dirac operator for the massless Wilson fermion. According 
to Ref. [^, locality of the overlap-Dirac operator can be studied by examining the locality of 
eigenvectors of Z/^. defined by H^(j)i{x) = Aw,,-0,(^). Since both the near-zero mode density and 
the corresponding eigenvectors are known to depend on the gauge action and j8, in this work we 
explore these two properties for several gauge actions. 

We consider three gauge actions, the standard plaquette (Plq), Iwasaki (RG), and a modified 
gauge action motivated by the admissibility condition (Adm). The Adm action [|3|, ^ ^] is given by 

l-ReTrPuvW/3 
3lruv : . , , when 1 - ReTrPu v/3 < e 



SAdn 



Otherwise 



where e is a parameter to control the possible maximum value of KeTrP^ vix) /3. With e < 
1/(6(2 + \/2)), the locality of the overlap-Dirac operator is guaranteed [^. This action reduces 
to the standard plaquette action when 1 /£ — > 0. In addition, for each of these gauge actions we 
introduce two-flavors of extra- Wilson fermions \j/ and ghosts ;^ as 

Sen = Y,¥{x)Dv,{-mQ)Y{x) + Y,x\x) (Dw(-mo) + /M/sT^) Z W, (1-3) 



where acts on the flavor index. Since integrating out these extra fields results in det \H^{—mo) 
/ {H^{—mo) + the appearance of the modes with |Aw,i| < /X are suppressed with this weight. 
jj. =0 corresponds to the standard quenched approximation. The same value of niQ is taken for both 
the overlap kernel and extra fields. 

The lattice size is fixed to 16^ x 32. For the Adm action, l/e = 2/3 is used throughout this 
study. Three values of jj. (0.0, 0.2, and 0.4) are examined, and more than 10,000 trajectories are 
accumulated for each gauge action as shown in Tab. |l]. The lattice spacing is set by ro=0.49 fm, 
and is about 0.125 fm for all lattices, unless otherwise stated. 



2. Spectral density 

Three panels in Fig. [l] show the distributions of the near-zero modes for three gauge actions, 
where all eigenvalues are plotted in the lattice unit. Different symbols represent different values 
of /I. We find that the extra fermions and ghosts suppress the low-lying modes as expected. The 
suppression is most effective with jj. =0.2. Figure ^ shows a comparison of the spectral density 
p(Aw) for the three gauge action with /i = 0.2. We observe that the RG action yields the smallest 
near-zero mode density. From these observations made in the quenched approximation, we decide 
to employ the Iwasaki RG gauge action with jj. = 0.2 in the dynamical overlap simulations. 
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action 


P 




l/e 


#of trj. 


Plq 


5.83 


0.0 




20,000 




5.70 


0.2 




11,600 




5.45 


0.4 




11,600 


RG 


2.43 


0.0 




20,000 




2.37 


0.2 




21,600 




2.27 


0.4 




20,000 


Adm 


2.33 


0.0 


2/3 


20,000 




2.23 


0.2 


2/3 


20,000 




2.06 


0.4 


2/3 


14,800 



Table 1: Simulation parameters. 
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Figure 2: Spectral densities for three quenched 
gauge actions, Plq (blue), Adm (black) and RG 
(red) with jx = 0.2. 
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Figure 1: Distribution of near-zero modes for three 
different gauge actions, Plq, RG and Adm from top 
to bottom, and with three values of /i=0.0, 0.2, 0.4 
from left to right. 
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Figure 3: Comparison of the near-zero mode den- 
sities on the quenched (black), slightly fine (blue) 
and coarse (red) dynamical configurations. 

We made the similar study on dynamical 
configurations, and compare the results with that 
on the quenched results with ;Li=0.2 in Fig. ||. 
The results from two dynamical ensembles with 
j8 = 2.30 and 2.35 are shown, their lattice spac- 
ings are a=0. 12 fm and 0. 1 1 fm, respectively, 
and the sea quark mass is about ms/4 in both 
cases. While the density increases for the dy- 
namical configurations, it is clear that the near- 
zero mode suppression works well even after 
incorporating the overlap sea quarks. More de- 
tails on our dynamical configurations are found 
in Ref. [p. 

3. Mobility edge 

The overlap-Dirac operator is proved to be 
exponentially local, if there is no low-lying mode 
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below some threshold [Q]. In practice, the density of the near-zero mode is non-zero, unless we 
introduce the determinant factor such as det//^. Even in this case, the overlap-Dirac operator is 
exponentially local if the near-zero modes themselves are exponentially localized Golterman, 
Shamir, and Svetitsky argued that the magnitude of the overlap operator behaves as 



\Doy{x,y)\ ~ Ap(A)exp 



-y\\ 



2l,{l)J 



+ 0(l)-exp(-A,|^-3;|), 



(3.1) 



where A denotes a near-maximum eigenvalue of low-lying localized modes (somewhat ambigu- 
ous), //(A) a localization length at A. The parameter Ac stands for the mobility edge, which sep- 
arates localized and extended modes. The first term is derived from the dependence of p ( A ) and 
Z/(A) on A, while the second term is a conjecture motivated by numerical experiences In most 
cases, it is known that the second term dominates the first. 

Mobility edge is determined from the eigenvectors as follows. We first define p,(x) and /,(r) 

as 



Pi{x) = 0/(x)0,-(x), Pi(xo) =max{p;(x)}, 

X 

fi{r) = {Pi{x)\r=\x-XQ\], 



(3.2) 
(3.3) 



where p, (x) is an average of p,(x) over the lattice points which have the same distance r from xq. 
The eigenvalues |Aw,i| are binned with a certain bin size, and fj{r) is averaged over the modes 
within a bin. The localization length at each bin, is then obtained at large r by fitting to 



=exp 



^/(Aw.i; 



(3.4) 



Mobility edge, A^, is defined by Aw,i at which //(Aw.i) diverges. Figure ^ shows an example of fj{r) 
obtained on a single configuration. We can clearly see that the decay rate becomes smaller for 
larger eigenvalues. 

The determination of A^ is performed with a fixed value of mo=l.6. In addition to the ensem- 
bles used in the study of spectral density, the determination is performed on three coarse lattices, 
which are generated with the Plq action with j3 =5.70 and 5.40 and the RG action with j3=2.43, to 



see the dependence on lattice spacing. //(Aw,;) extracted from the fit using eq. (3.4) is plotted as a 
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Figure 4: Example of /, (r). Different symbols represent different bin in Aw, as denoted in the plots. Data 
for the dynamical configuration at j3=2.30, OTsea=0.025. 
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Figure 5: Inverse localization length 1 /li{Xv,) as a function of \X, 



function of Aw,, in Fig. ^ for various ensembles. In the two plots at the top of Fig. ^, the results at 
the same lattice spacing, a=0. 125 fm, but with different jj. are plotted for the Plq (left) and the RG 
(right). Ac is found to be 0.24-0.27 for the plaquette action and 0.34-0.36 for the RG action in the 
lattice unit. Dependence of Ac on jj. turns out to be weak in the range of 0.0 < /i < 0.4. 

From the two plots in the middle showing the results with /i = for two or three different p , 
we find that Ac in the lattice unit decreases as jS decreases. The data at j8 = 5.40 for the plaquette 
action (circle in the left-middle panel) shows l///(0) ~ 0. Ref. suggests to use "l///(0) = 0" as 
the definition of the Aoki phase. 

In the left-bottom, the results with the plaquette and RG actions are compared at the same 
lattice spacing. Clearly the RG action gives lager Ac than that of the Plq action. Finally, the 
right-bottom figure shows a comparison of the results from quenched and dynamical runs with two 
different sea quark masses. Here we do not see any significant difference between the quenched 
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Figure 6: r dependence of h{r) for four different ensembles. 

and dynamical runs nor clear dependence on sea quark mass within dynamical runs. Ac turns out 
to be 500-600 MeV in our dynamical lattices. 

We also measure 1 /Zov from the r dependence of h{r) defined by 



h{r) =max< \\Do^iix,y)5{y - xo) 



(3.5) 



h{r) is measured on the same ensembles, and some of the results obtained with the RG action are 
plotted in Fig. ^, where numerical data are in the lattice unit. We see that all results except for the 
one with j8=2.20 coincide with each other. j3=2.20 corresponds to a ^ 0.2 fm, while others are 
(3=0.120-0.125 fm. By fitting this to the same form as eq. (3.4), we obtain 1 /Zov- For our dynamical 
lattice at j8=2.30, 1 //qv is estimated to be about 800 MeV. 

In Fig. ^ Ac and l/lov obtained with the RG action are compared. While both l//ov and 
Ac show the similar dependence on lattice spacing, the difference is sizable. Within the range of 
the lattice spacing we have studied it turns out that l/lov > Ac. Although Ac ~500-600 MeV in 
our dynamical lattice is not much larger than Aqcd> 1/^ov ~ 800 MeV is probably acceptable for 
simulations of QCD. In Ref. [^, l/lov has been studied at several different lattice spacings in the 
range of a ~0. 13-0.20 fm within the quenched approximation. They reported acceptably large 
values for 1 / /qv at all lattices, which appears to be consistent with our observation. 



•a 

•y 0.6 

-4—1 

^"0.4 

o 





RG, m^=1.6 



□ 1//^^,, (1=0.0 

o \n^^, (1=0.2 

A 1// , , Dynamical 



■ A.^ , (1=0.0 

A A. , Dynamical 



9 







0.1 



0.2 



Figure 7: a dependence of 1 //qv and A^ . 
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4. Summary 

In this work, near-zero mode density for various gauge actions are examined from the view- 
point of cost reduction in dynamical overlap simulation, and the RG Iwasaki action with extra 
fermions and ghosts exhibits smallest near-zero mode density among the actions studied. We also 
determine the mobility edge and the localization range of the overlap-Dirac operator. While Ac and 
I//0V seem to agree with each other qualitatively, the difference is sizable. 1 //ov tells us the locaUty 
of an given overlap-Dirac operator and is presumably all we need to know from the practical point 
of view, but the precise relationship between l//ov and Xc should be understood. 
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